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Abstract 

We explore the behavior of collective nuclear excitations under a multi-parameter 
deformation of the Hamiltonian. The Hamiltonian matrix elements have the form 
P{\Hij\) oc 1/ ^J\Hij\eyi\){—\Hij\/V), with a parametric correlation of the type 
\og{H{x)H{y)) oc —\x — y\. The studies are done in both the regular and chaotic 
regimes of the Hamiltonian. Model independent predictions for a wide variety of 
correlation functions and distributions which depend on wavefunctions and energies 
are found from parametric random matrix theory and are compared to the nuclear 
excitations. We find that our universal predictions are observed in the nuclear 
states. Being a multi-parameter theory, we consider general paths in parameter 
space and find that universality can be effected by the topology of the parameter 
space. Specifically, Berry's phase can modify short distance correlations, breaking 
certain universal predictions. 
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1 Introduction 



The statistics of nuclear excitations has been explored from the shell model to collec- 
tive models, with studies ranging from the relation of observed quantum fluctuations to 
those in random matrix models, to the connection with chaos using classical limits of 
the Hamiltonian[0-0. The agreement of various spectral properties with random ma- 
trix predictions has shown that certain simplifying assumptions can be made concerning 
fluctuations in nuclei. Once random matrix theory can be justified, certain results follow 
immediately. These studies of chaos in nuclei stem from attempts to extract a simpli- 
fied behavior from the complexity of nuclear excitations. In this respect, random matrix 
theory has provided invaluable assistance in developing simple methods to compute com- 
plex behaviors. In the past, aside from the studies of constant random matrices and 
the relation to chaos, these models have been given a parameter dependence to model 
correlations in various nuclear systems, from heavy ion collisions!^, high spin physics]^ 
to large amplitude collective motion|p. Recently it has been shown that Hamiltonians 
which have a parametric dependence can exhibit universal behavior]^, that is, there can 
exist model independent quantities in a given theory, providing the Hamiltonian has cer- 
tain random matrix properties. In this article we study a wide class of observables and 
develop universal predictions. We further show that parametric deformations of nuclear 
Hamiltonians can be readily modeled by a simple translationally invariant parametric 
random matrix theory, even though the Hamiltonian does not apriori look like a random 
matrix. We further justify the use of parametric random matrix theory for collective 
nuclear excitations. 



2 Collective Nuclear States 

We have chosen to model collective nuclear excitations in the framework of the Interacting 
Boson Model (IBM)||TD| for several reasons. One of our main objectives is to explore 
and categorize types of model independent predictions that exist in parametric quantum 
theories which exhibit classical chaos. The IBM is ideally suited for this since the classical 



limit has been extensively studied in recent years using coherent states [|ri|], and the 
complete chaotic behavior is now known for every value of the parameters [0. Hence 
we can easily choose parametric variations in regions of strong or weak chaos, or in 
regular regimes of the parameter space. An additional advantage is that we can solve the 
quantum problem exactly. One might argue that collective states form only a subset of 
the real spectrum as the excitation energy increases, so that the use of the IBM is not 
reasonable. This is not crucial, however, since the IBM provides a solvable theory with 
known spectral properties, which can be compared to those of the Gaussian Orthogonal 
Ensemble (GOE) throughout its parameter range. Certainly a more realistic description of 
the spectrum would embody the same features. For example, when broken pair states are 
added to the IBM model space, the spectrum becomes more GOE, as the interactions in 
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the Hamiltonian become more complicated [0]. This is certainly the case as one attempts 



to construct more realistic Hamiltonians. And as we are showing how model independent 
quantities emerge, the model we use is really not so important. Hence we use a simple 
form of the IBM Hamiltonian, known as the consistent-Q form: 

H = Eo + cifid + C2Q^ ■ + C3L ■ L, (1) 

where 

nd = S-d, = v^[dt X , Qx = d\s + s'^d^ + x[d^ ^d]f. (2) 

The parameters Cj are defined by ci = 1]/A and C2 = (1 — r])/4Ni,, where A^;, is the number 
of bosons. Since the Hamiltonian is diagonalized in a basis of fixed angular momentum 
L, the constant C3 does not play any role, and is hence omitted. Except when stated 
otherwise, we will use Nh = 25, which will give optimal statistics for the quantities we 
consider. The resulting dimensions for J'^ = 0''", 2+, 4"*", 10"'" states are 65,117,165,211. In 
this parameterization, one has the following limits: {a) rj = 1 corresponds to vibrational 
or U{5) nuclei, (h) t] = and x = — \/7/2 corresponds to rotational or SU{3) nuclei, and 
(c) = X = describes 7— soft or 0(6) nuclei. 

The interpretation of the Hamiltonian in terms of shape variables (3 and 7 is possible 
using coherent states, in the large limit of H. The energy surfaces for the Hamiltonian 



in Eq. (1) is p 



S{f3, 7; V, X) = f^-^ + - ^)(^ - 1) + cos 37^1 - ^(1 - r^)^. (3) 

For a particular value of t] and x, the energy E can be minimized to determine the 
quantities /3 and 7. /5 and 7 in turn define a deformed nuclear mean field. This can be 
made explicit by re-expressing the Hamiltonian in terms of excitations in a deformed mean 
field using boson condensate techniques 0]. This allows the interpretation of correlations 
of observables at different values of ri and x terms of the shape variable /3 and 7. 
Correlations in observables at different values of the parameters are then precisely the 
correlations between properties of the nucleus in the presence of different mean field 
configurations. We will consider the behavior of the properties of the Hamiltonian under 
very general parametric deformation z = z{ri,x)- For paths which lie entirely within 
the chaotic regime of the parameter space, the universal predictions we explore are path 
independent (up to effects due to Berry's phase which we explore in Sec. 5); correlations 
in a nucleus changing from rotational to vibrational or vibrational to 7— soft are the same 
when properly interpreted. 



2.1 Distributions and Correlations of Nuclear Matrix Elements 

One of the results presented in this article is that parametric nuclear Hamiltonians 
can be modelled by correlated, parametric gaussian random matrices. Recall that a 
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gaussian random matrix has a distribution of matrix elements of the gaussian form 
P{Hij) OC exp{—Hfj/2'y{l + Sij)), where 7 a constant related to the level density. To 
implement random matrix theory does not imply that the actual nuclear Hamiltonian (1) 
have gaussian matrix elements. We note that the distributions of matrix elements of the 
interacting boson model Hamiltonian are not gaussian. At any given value of {'r],x)y we 
find the distribution of matrix elements obeys roughly [|15 



P,U\H^^) OC ^=e-l^-^l/^ (4) 



where the strength V depends on whether one is in a chaotic or regular regime. Typical 
results are shown in Fig. 1 for both regular (crosses) and chaotic (boxes) choices of 
the parameters, together with the behavior (4) (solid curves). In the chaotic parameter 
regimes of the model, V is of order unity, while in the regular regions, it is much smaller. 
But both regular and chaotic regimes display the same functional form of the distribution, 
suggesting that the functional form is due to the structure of the Hamiltonian, rather 
than to the presence of chaos. Similar distribution functions have been seen in parity 



non-conservation studies of the compound nucleus 115 



Another quantity of interest is the autocorrelation function Fj^^ computed from the 
IBM Hamiltonian: 

F,U^-z') ^ {H{z)H{z')) = /^-^^^l-^^J2H^^{z)H,,{z')^ . (5) 

The brackets (■ ■ ■)z-z' represent the averaging over a trajectory in parameter space z{ri, x), 
which remains either in a chaotic or in a regular region, keeping the difference z — z' fixed. 
The results for the short distance behavior of the measured function Fj^^ are shown in 
Fig. 2, and are found to behave generically as 

i^i6m(^)~e-^l^l~l-7kl + --- (6) 

in both regular and chaotic regions. Again, the measured value of Fibm is not a good 
measure of the underlying chaos. If the averaging in (5) is restricted to a submatrix 
Ni < i,j < N2, there is no difference in the function Fif,m- One observation is that 
the apparent decorrelation, seen in the slow decay of Fif,m, is misleading. The actual 
decorrelation is much more rapid (when model specific dependencies are removed), as we 
will see below when we compute properties of wavefunctions and eigenvalues. Hence using 
Eq. (5) as actual input into a random matrix theory (e.g. into Eq. (16) below) is not 
meaningful. 



2.2 Unfolded Parametric Energies 

To study the statistical fluctuations of the nuclear energy levels, En{z), where z = z{ri, x) 
is a general path in configuration space, we must separate out the average behavior of 
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the energies which cannot be described by random matrix theory. This is done with the 
staircase function |Q 

N{E-z) = i:ie{E ~ H) = Y.e[E - En{z)] (7) 

n 

This function is computed along various paths in parameter space. On each path, 100- 
200 values of z are taken, and the eigenvalues {En{z)} determined. From this data, a 
polynomial fit is made to the staircase function using: 

k 6 

N{E■z) = Y.T.C^,z'E^■ (8) 

In the chaotic regions of parameter space. A; = 2 is sufficient, whereas higher values are 
needed is less chaotic regions. Once the coefficients Cij are determined, the unfolded 
energies are defined by 

Eniz)=N{En;z) (9) 

which results in a spectrum with a constant average level spacing. The degree of chaos 
in the energies can be measured through the Brody distribution of the level spacings 

S = En{z) - Er,.l{z)M 



P(s) = As"e-"^"^\ A={l + uj)a, " = ^ (^^) • (10) 

When u = 1, the distribution of level spacings is GOE, while for a; = 0, it is Poisson. 
In Fig. 3, we show the original (top) and unfolded (middle) parametric energies of the 
Hamiltonian for 77 = and x shown. The degree of 'chaos', measured in terms of the 
Brody parameter u, is shown in the bottom of the figure for this particular path. 

For purposes of contrast, the parametric levels (not unfolded) in two regular regions 
are shown in Fig. 4. The right figure corresponds to the (regular) transition from f/(5) 
to 0(6), while the left corresponds to a path through the 'valley of regularity' recently 
studied |]12| (the two kinks in the parametric energies are artificial and only refiect the fact 



that the parameter path took a turn). In the figure, the average Brody parameters ZJ are 
0.25 and 0.23, respectively, with fiuctuations up to = 1. In the following studies, we 
will consider regular and chaotic regions of the parameter space. As chaos is a classical 
notion and we study quantum statistics here, let us be precise. The chaotic regions are 
general paths in parameter space z{r], x) which stay in the areas of chaos in the classical 



limit of the Hamiltonian, as discussed in Ref. []T2|. The chaotic paths studied here are 
largely: 

i V = 0, XG [-0.9,-0.4] 

4v,x) = \ v = o.i, [-0.9,-0.4] (11) 

[ X = -0.66, r] G [0.0,0.5] 
Similarly, the regular regions we study are: 

^(^^.ix = -^V V^[0,l] 0(6) -[/(5) 

[0.5, 1.0] ^ ' 
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Also included in the regular region is the path through the valley of regularity from SU{3) 
to f/(5), as shown in the left half of Fig. 4. The first path in Eq. ( ]T^ corresponds to the 
right half of Fig. 4, is not the direct path between U{5) and 0(6) which would be x = 0, 
r] e [0, 1], and is entirely integrable. Rather we have chosen one that passes through a 
weakly chaotic regime. Hence, results in the following sections referred to as 'regular' 
are quantities which are averaged in the highest co regions of these regular areas, such as 
z e [0.5,0.8] in Fig. 4. 



3 Correlated Random Matrix Ensembles 

In order to see what types of model independent quantities emerge from the IBM, we 
must construct a random matrix model which has an equivalent parametric dependence. 
There is no unique method to realize such an ensemble, for instance one might take: 

H'{z) = H^ + zH2 (13) 
H"{z) = smzHi + coszH2 (14) 



H"'{z) = dyf{z-y)V{y), V{z)V{y) = 6{z - y). (15) 



Here Hi, H2 are constant N x N gaussian random matrices and Vij{y) is gaussian white 
noise for each i,j and y. Each of these is a viable random matrix theory, however the 
stochastic integral of Eq. (|15D, introduced by Wilkinsonfl^, provides a more general 
framework and includes a broader class of processes[0]. One additional difference be- 
tween H' and H", H'" is that the former is not a translationally invariant theory. While 
translational invariance is not important to the results we derive here, its presence sim- 
plifies our constructions. We will only focus on constructions of the type (14)-(15) here. 

The gaussian random matrix Hamiltonians are characterized by their first and second 
cumulants: 



H,,{z) = 



H,,iz)H,iiz') = ^F{z~z')g%, (16) 



where gfj^l^ = SikSji + SuSjk, gfjja'' = '^^u^jk, and a determines the average level spacing 
A through the relation a/A = \/2N /tt. Here z/ = 1 (GOE) corresponds to real-symmetric 
matrices, or equivalently, to a system with time reversal symmetry, and v = 2 (GUE) to 
complex hermitian matrices, or broken time reversal symmetry. From the definitions of 
H{z), it is clear that H{z) is GOE (z/ = 1) or GUE (z/ = 2) for any z. 

In contrast to previous studies of chaos in nuclei, which deal with constant random 
matrices, we can introduce a measure for the parametric ensemble rather easily ||18|| : 



P [H{z)\ oc exp j dzdz'Ti [H{z)K{z - z')H{z')]j , (17) 
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where the measure D [H{zy\ = Hz dH{z) is a product over the continuous variable z of the 
corresponding gaussian ensemble measure dH{z). Here Kij{z) can be viewed in general 
as a banded matrix of bandwidth a, connecting states i and j with |« — j| < cr. As we do 
not consider banded parametric matrices here, we will take Kij{z) = 6ijK{z), resulting in 
the measure: 

P [H{z)] oc exp J dzdz'K{z - z')Tt [H{z)H{z')] 

These gaussian integrals are easily done to establish Eq. (|16|), providing F is the inverse 
of K, 

' dz'K{z - z')F{z' - z") = 6{z - z"). (19) 



The stochastic integral (|T5D provides a direct method for constructing H{z) with a desired 
F{z). That is, we can choose / to satisfy 

F{z -y)= j dxfiz - x)f{y - x), (20) 

then H is constructed as in ([T5|), and the desired covariance (pIBf ) is automatically satisfied. 
It is important to realize that the properties of the random matrix theory here are distinct 
from the observed properties of the Hamiltonian, both in the measured distributions 
of matrix elements Pnym{Hij) and their autocorrelation Fibm{z). The random matrix 
distribution of matrix elements is gaussian, and F{z) is different from Fibm, as we discuss 
below. In particular, F{z) cannot be exponential as measured. Model independent results 
can be obtained from our random matrix constructions by a proper scaling of parameters. 
We discuss two approaches here to this scaling. The first is a general procedure based 



on Ref. [18], while the second is a more heuristic argument based on the Fokker-Planck 



equation and the original work of Dyson, and has been pursued in recent works|T^, 20, 21 . 



3.1 Scaling from Anomalous Diffusion 

It was shown recently that universal (model independent) predictions can be obtained 
from the above translationally invariant random matrix theories if one introduces a proper 
scaling! 18]. A general approach to do so is to view the parametric dependence of the 
energies as a diffusion process. Consider first the short distance behavior of the function 

F{zy. 

F{z) ^ 1 - c^\z\'' + ■ ■ ■ . (21) 
From perturbation theory, one can see that 



|2 



I r 77- 

5E^ix) = E^{x') - E^ix) = 5H„,^ + J2 ^ -E + - ^^^^ 

{^nix')\^n{x))\' = 1-E7#^^ + -- (23) 
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By implementing the ensemble averages defined in (16), and following Dyson [ll9|, one 
easily finds that: 

4/Vc 

{5E^iz)y ^ ^Sz^ ^ DJz'^ , (24) 



and similarly 1 — | (\[^„(2;) | \l/„(z')) p oc 6z°'. One can then view the parametric energy 
levels En{z), such as those in Fig. 2 (middle), as evolving diffusively on short distance 
scales according to Eq. (|2^. This has been recently contrasted with the anomalous 
diffusion process of a particle in a chaotic or disordered medium, whose position obeys < 
i?^(t) >= Df^, which although the physics is distinct, the formal treatment is similar p!8|. 
For our random matrix model, the parameter z{r], x) plays the role of time. The diffusion 
constant contains both dimensional information (A^) and model dependent data (cq). 
Hence, by scaling the parameter z by the diffusion constant D, all model and dimension 
dependence is removed. This is done by defining a new scaled parameter 

z=[D^Y/-z=(^^)'^\. (25) 



For the case of a = 2, we have D2 = C(0), where C(0) is the scaling introduced in Ref. 
By computing observables in the scaled variable, one obtains model-independent predic- 
tions for desired quantities. Physically, a = 2 in Eq. (^) corresponds to a Hamiltonian 
with a smooth dependence on the parameter z, while a < 2 corresponds to a theory 
with fractal parametric dependence, such as a parameter range taken from a Brownian 
trajectory. The computation of Da is done through the definition (|2^): 



Da = -^^^ C{0) = D, = (26) 

In order to model the IBM Hamiltonian using parametric random matrix theory, we 
must choose a correlator F{z) in Eq. ([16|) with a = 2 short distance behavior (see Eq. 
(HI))- Otherwise the parametric dependence of the random matrix energies En{z) would 
not be smooth. Hence, if we attempt to incorporate nuclear properties into the random 
matrix model by substituting the model specific, computed Fihm{z) in Eq. (6) into Eq. 



(p!6|), we would end up with an energy spectrum En{z) characterized by a = 1, resulting 
in non-smooth, brownian motion type paths for each En{z). For the IBM, a = 2 is the 
proper result. Details can be found in Ref. p 



3.2 Scaling and the Fokker-Planck Equation 



A more heuristic argument can also be made with the Fokker-Planck equation. Fokker- 
Planck methods, first introduced by Dyson|l^, have been recently discussed with para- 
metric correlations in mindp0|-||21||. Consider first Dyson's Brownian motion model for 
random matrices: 



Hij = --fHij + fij{t) 



(27) 
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The random force is white noise: 



Ut)fm = Tgl;l,6it-t'). (28) 

and 7 is a friction coefficient. The equihbrium solution can be found in the long time 
limit, by direct integration: 

Hij{t)= f e-'\'-^\h,{T)dT (29) 



which is same type of stochastic formulation as in Eq. (pISl). It follows then from Eqs. 
(0,0,^) that: 

{H,,mUt')) = ^^?£Le-^l*-*'l (30) 

This process can be formulated as well in terms of the Fokker-Planck equation for the 
distribution P{H,t): 

BP 3 1 f)^ P 

Since we are interested in a stationary process P{H,t) = P{H) for H{t), we can choose 
the initial distribution to be the equilibrium result P{H) oc exp[— z/Trif^/2a^]. The 
equilibrium solution is a solution of the Fokker-Planck equation providing the fluctuation 
dissipation theorem is satisfied: 

g^],T/2^ = IH-]^ = ^g^], (32) 

or a^/u = r/7, where the last equality follows from our original construction in Eq. ([161) . 
Following Dyson's argument that the eigenvalues of H{t) behave as a diffusive coulomb 
gas, equilibrating (at the microscopic scale) on a timescale of t oc l/i^jN), we see that in 
order to obtain A^— independent correlations, we must have 7 oc or 7 = 'j'/N. By 
equating the second cumulants of H(t) for and Langevin process (|30D and for our desired 
process ([16|), and implementing the fluctuation-dissipation theorem, we equate 

F{z-z')^exp(^-^\t-t'\^ . (33) 

Expanding to leading order in the large limit, we associate: 

1 -cl^r ~ 1 - — Itl (34) 

As the Langevin result is A^— independent at the microscopic scale, A^ independent results 
will also result if we deflne a new quantity z = {c/NY^'^z. Up to a factor of order unity 
(4/7r^z/), this is precisely the same scaling we discussed earlier. Now to leading order, 

Tf'^ \ Z -z' \"- 
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There is now no explicit model dependence c, and the A^— dependence of correlation 
functions will be absent due to the explicit dependence on N, which was require to achieve 
the microscopic equilibrium condition discussed above. 



4 Observables 

In this section we will see that many properties of parametric Hamiltonians have well 
defined model independent structure. There are two classes of observables we study here, 
those related to the energies En{z) and those related to the instantaneous eigenfunctions 
I The relation of wavefunctions | \l/„(z)) at z to those at z' is given by the 

transformation matrix: 

Unm{z - Z') = {^>n{z) \ ^m{z)) (36) 

We will see below that correlation functions that depend on Unm{z) and En{z) have 
universal predictions which generally agree well with the results of the IBM in the chaotic 
regions. 

The universal predictions for an a = 2, GOE system are computed here with two 
different covariances. The first is the simple sum of two uncorrelated GOE matrices 

H'{z) = Hi cos z + H2 sin z, H'{z)H'{y) = cos{z - y) = F\z - y). (37) 

The correlator is periodic with F'{z) ^ 1 — + ■ ■ ■ , defining the scaling z ^ z = 
a/D^z = {■\/2N /'k)z for universal correlations. The second construction is in terms of a 
stochastic integral, where one integrates over a continuous range of uncorrelated random 
matrices V{ii)\ 

H"{z) = J dye-^'-y^''I^V{y), H"{z)H"{y) = e-^^'yfl^ = F"{z - y) (38) 

Here z = {^/N/tc)z. We have computed various correlation functions described below 
with = 50 — 300. Generally = 50 is sufficiently large. 



4.1 C(z) 



The slope-slope correlation function of the unfolded parametric energies En{z) 
asi 

\ 9^ 9^ I E,z 



is defined 



(39) 



where the averaging is over energy and parameter. In Fig. 5, the results for the IBM in 
the chaotic regions are shown for J'^ = 0"*", 2"'", 4+, 10"^. The averages are computed by 
averaging over the middle third of the spectrum and over the trajectory z. The two solid 
lines are the results of the random matrix simulations, Eqs. (|37| - |38D , for = 50 and 
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300. For comparison, a computation in the regular region of the IBM is shown indicating 
much slower decorrelation. As we will see in all computations here, the typical distance 
at which quantities decorrelate is z ~ 1, which corresponds to the average separation 
between level crossings when the energies En are plotted as a function of the scaled 
parameter z. For regular systems, and the apparent level crossings in Fig. 4, this is not 
the case, and decorrelation happens over a much longer scale. In general the agreement 
with the universal functions is quite good. The 0^ states have the poorest agreement, 
which is in part statistical, having the smallest dimension. 



4.2 N-Scaling of C(0) 

We have seen that the diffusion constant C(0) = D2 scales linearly with A^, the dimension 
of the space. This scaling can be tested in the IBM by modifying the boson number A*";,. 
For A";, = 10, 15, 20, 25 we have dimensions of J = 10+ states of A^ = 16, 56, 121, 211. In 
Fig. 6 we plot C(0) as a function of A^, and see that this scaling is observed. At low Nb 
(equivalently low A^) , the results are not as reliable due to statistics becoming increasingly 
poor, and the classical phase space becoming increasingly regular |12|. 



4.3 Curvature Distribution P(k) 

The distribution of curvatures of the parametric energies En{z), or equally En{z), have a 
predicted distribution in the chaotic regime, given byp2|: 



Cu Cp En 1 

^(^) = (1 + ^5(^/2+1) ' ^ = (40) 

Here is the normalization, v = 1(2) for GOE(GUE), and k is the scaled curvature of 
the parametric energy. This function is compared to our random matrix computation 
(Eq. (^)) in Fig. 7 (a). In Fig. 7 (b), the results for the chaotic region of the IBM (solid 
histogram) are seen to agree equally well. A similar calculation done in the regular region 
shows a much more strongly peaked function (dashed, and scaled by 1/5 vertically). This 
is expected since the regular regions have much fewer level crossings and are hence much 
flatter (see Fig. 4). 



4.4 Diagonal Wavefunction Decorrelations: Pn{\Unn{z)[ 



The adiabatic survival probability \ {^n{z) \ \l/„(0))p measures how rapidly wavefunctions 
decorrelate. This was shown to be universal recentlyflS, 23|, with a well defined Lorentzian 
shape: 

f 1 Y 



Pn{\Unn{z)\') = \{^n{z)\^n{yM 



1 + c\z — y\ 



(41) 
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where v = 1(2) corresponds to GOE(GUE) eigenstates, c is a constant, and a is given 
by the leading order behavior of F{z). For the case at hand, u = 1 and a = 2. In 
Fig. 8, comparisons are shown for selected spins, and in chaotic (symbols) and regular 
(dashed) regions. The regular regions indicate much longer correlations, while the GOE 
result provides the most rapid statistical decorrelation of states. The two solid lines are 
the random matrix predictions. There is no equivalent universal prediction for the non- 
chaotic regimes, and the dashed line is just a representative Pn of the IBM in the regular 
region. 



4.5 Off-Diagonal Wavefunction Decorrelation: Pk{\Unm{z)\ 



2^ 



Wilkinson and Walker pl| have used perturbation theory to derive an approximate ex- 
pression for the distribution of squared off-diagonal matrix elements, \{'^n{z) \ \E'm(0))p, 
in the limit of \z\ oo and k = \m — n\ ^ 1. They found that 

Here, the energies are not averaged over, /i^ ~ 1, and p is the average level density. 
In exploring the behavior of these quantities in the IBM, it is difficult to satisfy the 
validity conditions for Eq. (|4^), since we largely study states in the middle portion of 
the spectrum, and m ^ n is difficult to satisfy. If we use only two states separated by 
= In — m| ^ 1, and we do not average over energy, the statistics are very poor. In order 
to get good statistics, we have examined the equivalent distribution which is averaged 
over both coordinate z and energy, with k kept fixed. We define this distribution of 
off-diagonal matrix elements as (A; > 0): 

Pk{z) = (\Unm{z)\^) = (\{^n{z + Zo)\^m{Zo))\^)^ , , , (43) 



E,Zo,k=\n—m\ 

where the subscript indicates that it is averaged over the trajectory z{ri, x), as well as over 
energy, with the separation \n — m\ held fixed. We use the notation Pnm for the quantity 
which is not energy averaged, and Pk for the energy averaged result. In Fig. 9, we compare 
this function computed in random matrix theory (solid) and in the chaotic regime of the 
IBM (boxes). In general the agreement is quite good. The random matrix theory result 
was done using N = 50, and averaging over the middle third of the spectrum. Hence as 
k increases, the statistics get worse. The IBM results are for the = 10"^ states, with a 
dimension of 211, so that the statistics is better. In order to contrast our results with Eq. 
(^), we have taken a rescaled form, which is not entirely justified, as the scaling by fi is 
distinct from D2 = C(0). Nevertheless, we plot in Fig. 9, the following rescaled functions 
(whose regimes of validity are indicated): 

P'kiz) = A;»l,|2;|»l (44) 

P''(^) = Prfxcos(A;x)e-^[^~^"P(-I^IW, fc>l (45) 
Jo 
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The function P^. (dots in Fig. 9) is an approximation of P^' (dot-dashed in Fig. 9) in the 
hmit 1^1 ^ 1, which was derived in Ref. As the value of k increases, there is better 
agreement with the exact results form the random matrix simulation and the IBM. We 
observe that a better overall fit can be found with the function (dashed line if Fig. 9) 



which trivially converges to the Wilkinson- Walker result in its regime k ^ 1, but better 
describes the results for all k. The results for the regular region x = —0.661, t] G [0.5, 1.0], 
are given by the crosses. Once again, there is no universal result for the regular case. 
Further, as expected, all three analytic functions (42)-(44) agree in the large k limit. 



4.6 Diagonal Matrix Elements: Pj{Unn{z)) 

The previous results have been averages over various matrix elements. We now show that 
the actual distributions of matrix elements can also be predicted by universal functions. 
Consider, for example, the distribution of the matrix elements Unn{z) = {"^niz) \ \E'n(0)). 
These can be seen to be described by a universal function for each z. The distribution 
PiJJnniz)) is shown in Fig. 10 for = 10+ states at several values of z. At z = 0, 
the distribution is a 5— function centered at 1. As the separation of the wavefunctions 
increases, the centroid shifts from 1 to 0, and the overlap evolves to an asymptotic Porter- 
Thomas distribution (an approximate gaussian) at large z: 



p,iuu^))^{ . ;r: (47) 



6{Unn - 1) Z^O 

exp(-Ara„V2) z>l 

The large z limit is not universal, as can be seen by the explicit dimensional dependence 
of the result. The Porter-Thomas result lies beyond the range of universality which is 
largely z < 1, but the precise decorrelation of matrix elements within this range has the 
model independent form shown in Fig. 10. 

4.7 OfF-Diagonal Matrix Elements: Pz{Unm{z)) 

A similar result exists for the distribution of off-diagonal matrix elements Unm{z) = 
{^n{z) I \E'm(0)) . The distribution P~{Unmiz)) is shown in Fig. 11 for J'^ = 10+ states 
at several values of z. At z = 0, the distribution is also a 5— function, but now centered 
aX z = 0. As the separation of the wavefunctions increases, the overlaps evolve to an 
asymptotic Porter-Thomas distribution (an approximate gaussian) at large z: 

i:: (48) 
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As before, the large z limit is not universal, as can be seen by the explicit dimensional 
dependence of the result. The solid curve in Fig. 11 is the Porter-Thomas result. 

4.8 Correlations of Mean Fields 

Each value of (r/, x) corresponds to a deformed mean field characterized by 7) de- 
termined from the minimum of Eq. (3). Because wavefunctions decorrelate on order of 
z = \/D2Z ~ 1, the actual correlation length in terms of the parameters rj and x depends 
on spin, and is given by 2 ~ = 1/ ^D2{A(3, A7, J'^). To explore the spin dependence 
of the correlation length, we compute z^. in the IBM for J'^ = 0^, 2+, 4+, 10^, and find 
typical values of z^ = 0.16, 0.14, 0.11, 0.05. This has a roughly behavior 

~ 1 - 7^ (49) 

where 7 is a constant. How generic such an dependence might be in other nuclear models 
is unclear, but it does indicate how rapid states of different spin can decorrelate (n.b. 
the results can be corrected for the dependence of Zc on ^/N, but this does not account 
entirely for the behavior). The equivalent values of A/5 and A7, which correspond to 
statistically decorrelated configurations, depend rather strongly on the parameter region. 
For example, in our calculations we can obtain a range from A/3 = 0.01 to 1.3, for the 
same correlation length, depending on whether the shape is undergoing a rapid shape 
phase transition or not in the particular parameter regime. One can only conclude that 
near a shape phase transition, there can exist strong statistical correlations between very 
distinct nuclear shapes. 

4.9 More Complicated Operators 

It is clear that one can explore many classes of operators and establish the behavior of 
model-independent limits of those quantities. For instance in the study of the E2 decay 
of high-spin states, Aberg^ has introduced the matrix quantity 

T,, =1 (vl/,(j - 2) I ^,iJ)) |2 (E.iJ) - E,iJ - 2))^ (50) 

where the parameter J is the angular momentum, and is equivalent to z. This matrix can 
be explored as a function of the correlation length, and has different results in the chaotic 
limit, depending on the spin dependent scaling One can consider other operators as 
well, and we would like to point out that additional quantities can be constructed using 
our universal predictions here, together with the analysis of Ref. [^, which discusses 
how to compute arbitrary correlation functions. 
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5 Multi-parameter Correlations: Topological Effects 
and Berry's Phase 



While formal studies of parametric correlations have been limited largely to single param- 
eter systems (see Ref. for some exceptions), nuclear deformation is usually described 



in terms of two or more shape parameters. When two or more parameters are involved, 
one finds that short distance correlations can be modified by topological effects, due to 
Berry's phase. That is, the correlation between quantities at /5, 7 and 7' depend on 
the path used to connect these points. Generally, for correlation functions which are sen- 
sitive to phase information, we will show that interference terms can strongly modify the 
expected results. We explore the basic ideas here in the case of two parameters. 

When a wavefunction undergoes parametric evolution on a closed circuit C, it is well 
known that the wavefunction can pick up a topological phase: 

m^^^)^e'^ic)^n{z), (51) 

where C represents a loop in parameter space starting and ending at z. For real symmetric 
matrices, such as our GOE ensemble, 7(C) is only or vr (mod 27r)p4[. Hence 



^^{z + C) = ±^r,{z) (52) 

where the sign depends upon the particular eigenstate and the path, and z + C represents 
the same point z after following the closed loop C. Of course, one does not have to follow 
a closed loop. A similar effect exists if one follows two distinct paths from z to z' . Then 
phase differences result in interference. Whether or not paths are in a chaotic or regular 
regime does not change the flavor of the argument, but in the chaotic regime, more states 
are likely to pick up a negative phase due to the many avoided level crossings [B^]. 



5.1 A 2-Parameter Random Matrix Model 

The simplest formulation of a two parameter correlated random matrix ensemble is 

H{X, Y) = Hi cos X + H2smX + H3 cos Y + sin Y (53) 

where the constant random matrices Hi are uncorrelated: HiHj = 6ij. It follows that 

H{X,Y)H{X',Y') = cos(X - X') + cos(F - F' ) . (54) 

Generalizations to arbitrary dimensions have been discussed by Wilkinson We can 
now consider parametric excursions in the (X, Y) plane, specifically two paths which 
connect (0, 0) to (AX, AY), one of shortest length, and the other a longer path enclosing 
an area A. Because the wavefunctions acquire a Berry's phase around the closed loop, 
which can be ±1 for the GOE case, the area A enclosed can modify expected short distance 
behaviors. 



15 



5.2 Correlations in the /3 — j plane 



We can now explore some of the topological effects in our two parameter theory H{r], x)- 
Consider a rectangular loop C in parameter space which encloses an area A. In analogy to 
scaled parameter z, we define the scaled area of the loop a.s A = C{0)A ~ z^. Then an area 
A ~ 1 is a loop whose sides are approximately the decorrelation length of observables 
Such a loop stays within the universal regime for all values of the parameter. In Fig 
12, we plot the distribution of matrix elements PziJJnn) (see Eq. (45)) for such a loop 
Starting from the top of the figure, we have z = 0, and the distribution is a delta function 
As z increases, the distribution spreads in accordance with universal predictions (cf. Fig 
10). At the farthest point of the loop, the distribution is given by the middle figure. As 
the trajectory returns to the initial point, approximately half the eigenfunctions develop a 
negative topological phase, and at the final point, which is precisely the initial point, the 
distribution is equally split. All of the results in Fig. 12 are within the universal regime, 
but one can see that topological effects can destroy the expected behavior discussed in 



Fig. 10. For smaller loops, the effect is smaller. The approximate behavior is[p5 




/ = <i ? (55) 



Here / is the fraction of the total states which split to —1, and A is the enclosed, scaled 
area. The fraction increases linearly with the area. Because saturation occurs near A = 
1 = C(0)AxA?7, and C(0) oc A^, the size of the loop needed to see the maximal effect 
decreases like 1/A^: A^Ar/ ~ 1/A^. Hence we see that universal predictions can be 
modified in multi-parameter theories due to the topology of the parameter space. It is 
not sufficient to give only the metric distance in parameter space in order to provide 
all universal predictions. One must also consider the path taken to get to that metric 
separation. 

So there are several aspects here to consider. Berry's phase effects are independent 
of the underlying chaos along the chosen path, but depend more upon the nature of the 
parameter space enclosed by the path. One could imagine a loop in parameter space 
which is entirely regular, but encloses a chaotic regime. The expression for the fraction 
of total states / above assumes one is always in a chaotic regime, and the enclosed area 
is also chaotic. 



6 Conclusions 

In conclusion, we have explored the adiabatic behavior of collective nuclear excitations, 
and found that under the appropriate scaling of the parameter, correlation functions 
and distributions of matrix elements behave universally. Hence, if we wish to implement 
random matrix theory to study a complex nuclear situation, we specify immediately a 
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multitude of model independent results related to the wavefunctions and energies. The 
results here indicate that a new universality exists in nuclei, related to the 'deformation' 
of the nucleus, which is quite robust. As the random matrix predictions are generic, 
they should be present in other classes of nuclear states, generated from the shell model, 
or other models. While we have focused on matrix element distributions and certain 
correlation functions, it is clear that the scaling provides a general type of approach to 
compute arbitrary correlation functions. This also establishes that the use of random 
matrix theory with covariances of the type (16) is quite reasonable. 
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Figures 



Figure 1. Distribution of matrix elements for the IBM in the chaotic region (boxes, 
7] = 0,x = —0.7), and in the regular region (crosses, r] = 0.85, x = —0.661). Both 
distributions are compared to the distribution P{\H\) oc lyJ\H\ exp{—\H\/V), with values 
of y = 10 and V = 0.3, respectively. 

Figure 2. Short distance behavior of the measured IBM autocorrelation function {H{z)H{0)) 
in the chaotic region (solid) and in the regular region (dots). The linear behavior at small 
z suggests a = 1 in Eqs. (21) and (26), which is not consistent with the observed para- 
metric energies En{z). Hence Fibm cannot be used as physical input into the random 
matrix theory through Eq. (16). 

Figure 3. Instantaneous eigenstates of the Hamiltonian (1) for the parameter range rj = 
and X ^ shown. (Top) Original energies; (Middle) unfolded eneries; (bottom) Brody 
parameter along this path, indicating a rather chaotic regime, u = 1 being the GOE 
limit. 

Figure 4. Instantaneous eigenstates of the Hamiltonian (1) for two largely regular regions. 
(Left) A path from rotational {SU{3)) to vibrational (f/(5)) spectra through the regular 
region proposed recently [|T^ . The energies have been scaled by 1/2. (Right) The transition 
from vibrational (f/(5)) to 7— unstable (0(6)) choosing a path which is weakly chaotic. 
The average Brody parameters are ZJ = 0.23 and cJ = 0.25, respectively. 

Figure 5. Slope-slope autocorrelation function, C{z), for the parametric energies Eniz) in 
the chaotic (symbols) and regular (dots) regions. The solid lines are the random matrix 
predictions using (37) and (38) with = 300 and N = 50, respectively. The 0"*" states 
have the poorest statistics, due in part to the small dimension of the space (N=65). 

Figure 6. Scaling of D2 = C{0) with N, the dimension of the space. The boson number 
was varied as A*";, = 10, 15, 20, 25, resulting in dimensions A^ = 16, 56, 121, 211 for J'^ = lO"*" 
states. The anticipated scaling behavior, given by Eqs. (25)-(26) with a = 2, is linear, 
shown by the solid line. There are deviations at small boson number since the chaos is 
not as strong there, and the dimensions are small. 

Figure 7. (a) Analytic level curvature distribution P{k) (solid) compared to results of a 
GOE simulation (histogram), (b) Comparison of the analytic distribution to those for 2"*" 
states in the IBM (solid histogram). The dashed histogram corresponds to 2+ states in 
the regular regime, and has been scaled vertically by 1/5. 



Figure 8. Wavefunction decorrelation function Pn{z) = |(^„(z)|\I'„(0))p for selected 
states in the chaotic regimes indicated (symbols), and one for 2"^ states in a regular region 



(dashes). The solid curves are our universal predictions |jT8|. As there are no universal 
predictions in the regular regimes, the dashed curve is only representative. 

Figure 9. Distributions of off-diagonal matrix elements Pk{Unm{z)) (where k = \n — m\). 
The IBM results for chaotic J'" = 10"^ states (boxes) agree well with the random matrix 
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predictions (solid) using = 50, as well as with a simple analytic function (dashes). In 
addition, we show the asymptotic results Pkiz) (dots) and Pk{z) (dot-dashed) of Ref. 
which converge for low k to the exact results. For reference, a similar calculation in the 
regular region for 10+ states is included (crosses). 

Figure 10. Distributions of diagonal matrix elements P~{Unn{z)) at several values of z for 
lO"*" state (solid) and the random matrix predictions (dashed). The distribution shifts from 
a delta function centered at one at ^ = 0, to an asymptotic, non-universal Porter-Thomas 
distribution for ^ ^ 1 

Figure 11. Distributions of off-diagonal matrix elements P~{Unm{z)) at several values of z 
for lO"*" state (histogram). The distribution shifts from a delta function centered at zero 
at z = 0, to an asymptotic, non- universal Porter-Thomas distribution for z ^ 1 (solid). 

Figure 12. Effects of Berry's phase on universal distribution functions. The parameter z 
undergoes motion on a closed loop, starting and ending at ? = 0. Starting at the top, the 
distribution evolves and eventually bifurcates due to the presence of topological phases. 
Both results z = 0.5 are within the universal regime, but the lower figure shows that the 
path taken to get to the point of interest can be important. At the bottom, the fraction of 
matrix elements / that change to Unm = ~1 depends on the area enclosed. The minimum 
area needed to achieve the maximum fraction of 1/2 scales as 1/A^, and is hence rapidly 
realized in large dimensional systems. 
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